Calculation of Atomic Number States: a Bethe Ansatz Approach 
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I. INTRODUCTION 

The realization of Bose-Einstein condensation (BEC) 
in dilute gases has enabled the study and control of many- 
body systems. While most of the work has focused on the 
properties and excitations of the condensate, it has pro- 
vided a new path towards generation of atomic number 
(or Fock) states. These few body states with a definite 
number of atoms in the ground state are of great interest 
for quantum information where individual qubits can be 
addressed [H, @, Q , and could also be important for atom 
interferometry in order to reach the Heisenberg limit of 
detection. Experimentally, realization of Fock states re- 
quires a BEC confined in an optical box coupled with 
single-atom counting. The challenge is to obtain confine- 
ment in a trap that is comparable to optical lattices, but 
with only a single site [i| . Recent experimental work has 
demonstrated all the necessary steps towards this goal 
and they are now being incorporated into a single system 
In parallel, the theoretical analysis of this problem 
has focused on the conditions for optimum number state 
production [f| @] • These include the role of varying trap 
depth and size, either separately or in tandem. The in- 
teraction strength is a third control parameter that can 
be tuned with the the help of Feshbach resonances or 
by tuning the transverse confinement of the optical trap 
[j, Q and is considered here in more detail. 

It is clear that strong repulsion between atoms are 
desirable for the production of number states. The in- 
finitely strong interaction regime is the so-called Tonks- 
Girardeau regime where calculation has been made trivial 
thanks to the boson-fermion correspondence 0]. How- 
ever, this is only an unreachable limiting case. To make 
realistic predictions about what can be experimentally 
realized, we must consider the regime of relatively but 
not infinitely strong interactions. Previous studies of this 
regime were carried out using numerical methods which 
are both time-consuming and reveal little insight on how 
number states vary with interaction strength. In this 
paper we develop an approach that takes interaction as 
a parameter and by which we are able to chart number 



states in the parameter space with interaction as one of 
its dimensions [10| . 

Our analysis of the production of atomic number states 
in a ID optical box is based on the so-called Bethe ansatz 
approach. This method was first developed by Hans 
Bethe to solve the problem of a one-dimensional (ID) 
spin 1/2 Heisenberg ferromagnet [ll|. Since it's inven- 
tion, Bethe's method has found important applications 
in the study of interacting spin systems [H, 0, E3- 
It has also been applied to solve the problem of a 
ID bosoni c g as with repulsive (5-function interactions 
E EE EE El El HI HI • An outline for the structure of 
this article is as follows. In section UH we formulate the 
problem; in section ITTTI we present approximate solutions 
to our problem with Bethe ansatz approach; in section 
TV! we use Bethe ansatz solutions to analyze issues re- 
lated to number state production. 



II. FORMULATION OF THE PROBLEM 

The problem of many bosons with a (5-function interac- 
tion trapped in ID square well potential with finite well 
depth was studied using the Bethe ansatz in Ref [2l| . 
We find the problem of producing Fock states in the ul- 
tracold atom systems trapped in ID optical box bears 
similar characteristics. We treat the ID optical trap as 
a square well potential of length L and depth Vq. We 
write the interaction potential as ^c5(xi — Xj), where 
Xi and the positions of the interacting particles, H 

is Planck's constant and m is an atom's mass, and c is 
the interaction strength and has dimension of [1 /length] . 
According to [22j we have the following expression for c, 
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where a is the s-wave scattering length in 3-dimensional 
space, a± = \/2h 2 / 'mu; ~±, uj± is the transverse trapping 
frequency, and C ~ 1.4603 is an empirical constant num- 
ber. Since the interaction strength c depends on both 
scattering length and transverse trapping frequency u>±, 
tuning either of them will affect the interaction strength. 
The transverse trapping frequency may be controlled by 
optical box parameters [4| . Scattering length may be ad- 
justed by Feshbach resonance p|. To give a sense of order 
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of magnitude, for sodium atoms trapped in a ID optical 
box with transverse trapping frequency u>± = 2n x 150 
kHz and zero magnetic field, we have c = 16863.6 cm -1 . 
For 87 Rb atoms in a similar trap with zero magnetic field, 
we have c = 92391.6 cm -1 . 

To make our equations dimensionless, we use 1 /c as the 
length unit and h 2 c 2 /m as the energy unit. The square 
well potential is then 

v{x)=\ - k °y 2 > {x { < x °/ 2 ' (2) 

[ 0, otherwise, 

where fco and xq are dimensionless numbers. With these 
parameters, the well width is L = xq/c and well depth is 
V = h 2 c 2 k Q 2 /2m in cgs unit. 

The hamiltonian for the many-body system may be 
written as 

1 n „ 2 N N 

^-aE^+E^+E 5 ^-^)- ( 3 ) 

2—1 1 2—1 2J — 1 

i>j 

Our first main step is to solve the following eigenvalue 
problem 

%i>{ ~x) = Eip(~x), (4) 

where ~x is shorthand for xx, %2, • • • ,xn- We are primar- 
ily interested in bound states whose wavefunctions must 
be normalizable. As a minimum requirement, the wave- 
function of a bound state must satisfy lim^^-too ip(~x) = 
0. 



III. BETHE ANSATZ SOLUTIONS 

In section [TTl we defined a boundary value problem 
relevant to the production of number states. We will in 
this section apply Bethe ansatz and obtain solutions for 
it. 

As studied in previous literatures [H, E3, IHj], Bethe 
ansatz solution of this problem introduces a set of N 
as-yet- unknown wave numbers k — {kx,k2,- • ■ , &at}. 
In conjugate to these wave numbers, another set ~k = 
{ki, K2, ■ ■ ■ ,kn} is defined as 

«3 = \f k 0~ k 3 2 > ( 5 ) 

for j — 1,2, ••• ,N. The total energy of the Bethe 
ansatz state is E = — Y^j Kj 2 /2 = (k 2 — kg 2 ) /2. 
The eigenfunction (wavefunction) of Eq.(f5|) is piece- 
wise continuous in the TV-dimensional coordinate space 
{xx, %2, • ■ ■ , xn}- For simplicity, we consider three repre- 
sentative regions in the TV-dimensional coordinate space: 

Ri : —x /2 < xx < X2 < ■ ■ ■ < x^ < Xq/2, (6) 
i?2 : Xx < —xq/2 < x 2 < ■ ■ ■ < xn < xq/2, (7) 

-^3 : —Xq/2 < Xx < ■ ■ ■ < XN-1 < Xq/2 < Xpf. (8) 



Rx represents a region where all particles are trapped; 
^2,(3) represent a region where the 1st (iVth) particle 
tunnels into the left (right) barrier. In fact, each of these 
regions falls in a class consisting of AH regions that are 
related by coordinate permutation. For ease of reference, 
we name A the class of regions that can be obtained 
from i?i by mere coordinate permutations and study the 
wavefunctions in these regions at once. 

We denote the wavefunctions in a region of A as 
(j> T (~x), where r is the permutation operator that trans- 
forms R\ into this region, i.e., R T = tRx- This wave- 
function is the superposition of pure plane waves with 
"± signs" times permuted wave numbers, 

E E A{^a;T)e^^. (9) 

where q = {q,,^-" ^n} represents a possible combi- 
nation of N signs each of which is either + or — and 
C2 represents the group of such operations, G is the 
permutation group of N particles, and A(s, a; r) is the 
superposition amplitude. 

It is clear that the superposition amplitude A{<;, a; r) is 
a functional of the sign-flipping operator, the wave num- 
ber permutation operator, and the region permutation 
operator. By bosonic particle permutation symmetry, we 
establish the first set of equations among the superposi- 
tion amplitudes, 

A(s ) a;T)=A(s 1 ra;T), (10) 

where I is the identity element in the permutation group. 

The wavefunctions in region R2 have a more compli- 
cated form, 

2 (^) = E E B(<;,a)e {(Tli ^ xl e i ^f=^" at) ^, (11) 

? o-GG 

where (<^er k )j is the jth component wave number after 
the permutation operation a and the sign-flipping oper- 
ation <r, (a r K) 1 = \j k"^ — (cr k )i , which can be regarded 
as an extra operator on top of the permutation operator 
cr, and B{q,a) is the superposition amplitude. Similarly, 
the wavefunctions in region R3 may be written as, 

faCit) -EE CfaayEfe 1 (^)i^e-^^, 

(12) 

where (ct~k) n = \J k 2 — (cr k ) jy and C{q, cr) is the su- 
perposition amplitude. 

From Eq.((5]), it is clear that if, for any i, fcj > kg then 
there will be a corresponding pure imaginary Kj. From 
Eci. tTTTj) , Eq. (ri"2"|) , any pure imaginary Ki will cause the 
Bethe ansatz wavefunction unnormalizablc. Then from 
the normalizability requirement stated at the end of Sec- 
tion n we reason that a Bethe ansatz state is bound if and 
only if all of the wave numbers are real and smaller than 
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Icq. Since we are primarily interested in bound states, 
from now on we implicitly mean bound state when we 
say Bethe ansatz state, unless otherwise stated. 

Once we get the Bethe ansatz wavefunctions written 
down, the rest is straightforward. The main features of 
Eq. are the singular <5-function particle-particle in- 
teraction and the nonzero potential step at the edge of 
the square well. Bethe ansatz method elegantly treats 
both as boundary conditions. The boundary conditions 
at xi — Xj = for i, j — 1, 2, • • • ,N in regions of class A 
requires continuity of wavefunctions on the one hand, 



(13) 



and certain discontinuity in their first-derivatives on the 
other, 



dip dip 
dxi dxj 



dip dip 
dxi dxj 



The boundary conditions at Xi = ±xq/2, i = 1, 2, ■ • • , TV 
require continuity of both wavefunctions and their first- 
derivatives: 
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dip_ 
dxi 
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(15) 
(16) 



=x /2" 



Plugging Eq.©, ([HD, and (HU) into Eq.JTSD, d, 
and (|16p and including Eq. (fT0)) . we obtain the complete 
group of equations for our original problem. 

For the purpose of this article, it suffices to keep just 
the wave numbers and eliminate all other unknowns. Do- 
ing so yields the following secular equations for the N 
wave numbers, 



—Xokj — —nlj + 2 sin 1 I — 1 + 

N 

+ ^2 [ tan_1 ( k j + k i) + tan_1 ( k j - fez)] > ( 17 ) 

where j — 1, 2, • • • ,N, and / = I2, ■ ■ ■ , In} is a set 
of preselected integers. 

Now a retrospect on the three regions we selected in 
the coordinate space might make you wonder why we 
haven't included more regions (and more equations), pos- 
sibly with 2 or more particles lying outside the trap area 
simultaneously. In reality, this is possible, but only with 
small probability for deeply bound state. The effect on 
the number-state condition should be ignorable. Firstly, 
insofar as representative limiting cases (strong interac- 
tion limit and deep trap limit), Bethe ansatz solution 
agrees with known results (see Fig. [5h.). Secondly, the 
energy spacings between single-particle levels are big near 
the strong interaction regime where number state experi- 
ments take place most likely and therefore the probability 
for more than one particle to tunnel into the barrier is ig- 
norably small at all times. For this reason, we argue that 
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FIG. 1: (color online) Bethe ansatz states for sodium atoms. 
The total energy of 4-particle bound states are plotted against 
trap depth. The numbers at the beginning of each energy level 
are the quantum numbers of the bound state. We shifted the 
energy zero to the bottom of the trap for ease of presentation. 
Trap size L = 5 fim. Transverse trapping frequency u>± = 
2tt x 150 kHz. 



the Bethe ansatz-based approach is a sufficiently good 
approximation to our problem. 

There are a few notable aspects of Eq. (fT7)) . Firstly 
it is a transcendental equation and so does not have an- 
alytic solution. Secondly, we must pick a set of integers 
/ = {I±, I2, ■ ■ ■ , In} before we start the numerical com- 
putation. Apparently, not any set of integers would lead 
to a physically meaningful solution. Then a natural ques- 
tion to ask is what does. As argued in Appendix [A1 we 
find that Eq. p7[) yield valid solution if and only if the 
set of integers I are mutually distinct, somewhat similar 
to the theorem in Ref. [1 7]. A corollary of is that wave 
numbers thus obtained have mutually distinct absolute 
values. Because of this one-one correspondence, we use 
the set / as the quantum numbers for the corresponding 
Bethe ansatz state. 

Besides knowing what set of integers are valid quantum 
numbers, we need further identify the ground state and 
the first-excited state. As argued in Appendix [Bl we find 
that the ground state of an iV-boson system has the quan- 
tum number {1, 2, • • • , TV}; the first-excited state has the 
quantum number {1, 2, • • • , N - 1, N+ 1} (see Fig. [TJ). 

Without loss of generality, we re-organize the set of 
wave numbers such that < ki < ki < ■ ■ ■ < fcjv < feo- 
We define the jth single particle energy as Cj — —Kj 2 /2, 
for j = 1,2, ••• ,7V. Fig. [TJ gives total energies of the 
low-lying Bethe ansatz states. 



IV. NUMBER STATES 

We now apply the results of the previous sections to 
the problem of number state production. We can safely 
assume no atom with positive energy presents in space 
near the trapped area. In reality, if an atom acquires pos- 
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FIG. 2: (color online) Single particle energies of 4 sodium 
atoms in Bethe ansatz ground states. Trap size L = 5/wn. 
a. Dependence on interaction strengths (c). Trap depth Vo = 
kB x25nK, where ks is Boltzmann's constant. The dotdashed 
vertical line denotes the maximum interaction strength above 
which no Bethe ansatz state of 4-boson system exists. The 
other two vertical lines denote the interaction strengths of 
sodium (dotted) and 87 Rb (dashed) atoms at lo± — 2ty x 
150 kHz and zero magnetic field. Inset, the trap depth is 
lifted to ks x 40nK, at which condition all 4 atoms remain 
trapped to the Tonks limit, b. Dependence on trap depths 
(Vo). Transverse trapping frequency ui± — 2n x 150 kHz. 
Magnetic field is zero. The vertical line (dotdashed) denotes 
the minimum trap depth below which no bound state of 4- 
boson system exists. 




3.47 10.4 18.3 27.2 37.1 

Vo /(ft 2 /mL 2 ) 

FIG. 3: (color online) Ionization thresholds of sodium atoms 
with all parameters fixed except trap depth (we use h 2 /mL 2 
as energy unit for convenience). Only highest single particle 
energies of the Bethe ansatz JV-boson states are shown with 
N — 2 (circle), 3 (square), 4 (diamond), 5 (upright triangle), 
and 6 (invert triangle). Trap size L = 5 /im; Transverse trap- 
ping frequency ui± — 2ir x 150 kHz. The ionization thresholds 
(with the current numeric calculation step size) are marked 
with vertical lines (dot-dashed). 



itive energy, it would be quickly swept out of the vacuum 
chamber. Therefore, it is safe to assume that the states 
in the continuum spectrum are virtually unoccupied all 
the time. 

As dictated by Bethe ansatz, for some given trap pa- 
rameters (depth Vo, trap size L, scattering length a, and 
transverse trapping frequency ui±), N bosons can be con- 
tained in the trap if and only if there is an JV-boson 
Bethe ansatz state. The energy levels for an JV-boson 
system has been computed by Eq. (|17p numerically. As 



an example, Fig. [2] shows a 4-boson Bethe ansatz state 
that ceases to exist in certain regions of the parameter 
space. The Bethe ansatz state can only exist down to 
a certain trap depth in panel a and can only exist up 
to certain interaction strength in panel b, provided that 
all other parameters are held unchanged. We call the 
maximum number of particles that can be contained in 
the trap the trap capacity. The trap capacity puts an 
upper bound on the number states for a given point in 
the parameter space. The whole parameter space is thus 
partitioned into zones of certain trap capacities. We de- 
fine the boundaries of these partitions as the ionization 
threshold, because as we cross the boundaries from higher 
trap capacity side to lower capacity side adiabatically, the 
system state changes from stable to unstable and must 
release some particle(s). In Fig. |3]we show the ionization 
threshold for a system of 2, 3, 4, 5 and 6 bosons. 

Now that we have a clear upper limit, the trap capacity, 
on the number state that may be present for a given trap 
and other physical parameters. It remains questionable 
whether or not the trap capacity can be reached. The 
adiabatic laser culling technique developed in Ref. @ 
and the simulations in Tonks-Girardeau region made in 
Ref. [3] seem to suggest that it is possible to reach the 
trap capacity with the ultra-cold technique developed in 
Ref. H||. 

The starting point is a almost-pure Bose Einstein con- 
densate (BEC) that is optically trapped. Ignore ex- 
citation effects for now, it is useful to view the pro- 
cess from the angle of quantum optics and regard the 
state of the BEC as a coherent state (23[. A coherent 
state is essentially a superposition of Fock states with a 
Poisson-distribution in the boson numbers. As we adi- 
abatically change experimental parameters from parti- 
tions with higher trap capacity to lower capacity target- 
ing some number state, the Bethe ansatz solution puts a 
tighter and tighter restriction on the maximum number 
state. The system state thus undergoes two changes side 
by side: 1, more and more high-energy atoms are forced 
out; 2, more and more high-number Fock states are quan- 
tum mechanically 'projected' out of the system state (re- 
sulting in the so-called squeezed state). Each of these two 
changes has its distinctive effect on the system state: the 
first leads to smaller and smaller average particle num- 
ber N = (N) whereas the second leads to a reduction in 

the number uncertainty a 2 = (^N 2 — N J. Under opti- 
mal experimental conditions, the process continues until 
at some point, while the average number (N) > 0, the 
number uncertainty a « 0. A rigorous simulation of this 
would require calculating the value = as a function of 
time in a dynamic process and is certainly beyond the 
scope of this article. In Fig. [4] we show trap capacities 
and ionization thresholds as functions of trap depth and 
size. The interaction strength is implicit in the unit we 
adopted. 

There are several ways to tune the physical parameters 
to achieve the above goals. In previous references [1,0] , 
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FIG. 4: Map of number states and the ionization thresholds 
for sodium atoms in ID optical trap in adiabatic limit. Trans- 
verse trapping frequency lu± = 2tv x 150 kHz and zero mag- 
netic field are assumed. Note that c _1 and h 2 c 2 /m are used 
to make the axes dimensionless. a. Contour plot of number 
states as function of trap depth and size; b. and c. cross- 
sectional cut views at the indicated trap size (5pm) and depth 
(4nK) (solid lines) of a. The ticks on horizontal axes give the 
calculated ionization thresholds. 



only culling (changing trap depth), squeezing (changing 
trap size), or some combinations of the two are discussed. 
In certain circumstances, we propose that atom-atom in- 
teraction strength c be possibly tuned to supplement the 
production of number states. In view of the intrinsic lim- 
itations in tuning the trap parameters, it's possible that 
tuning of interaction strengths could play a key role in 
number state experiments. 

The path to a number state becomes clear now. By 
tuning the physical parameters of the ID optical trap adi- 
abatically, we force the ultracold atom sample through a 
series of quantum collapses until it eventually reaches the 
desired Fock state with some acceptable fidelity. Ideally, 
the course connecting the starting point and a targeted 
Fock state consists of a series of states (the Bethe ansatz 
states) with well-defined particle number. But in reality, 
there are always some elementary excitations, which is 
defined as any deviation from the ideal adiabatic course. 
Possible elementary excitations include occupations of 
excited Bethe ansatz state (of the same particle num- 
ber), earlier ionizations (loss of particles before reaching 
the Bethe ansatz ionization threshold) , and simultaneous 
ionizations of more than one particles. 

We now analyze the effects of excitations. Abrupt 
changes in the trapping potential tend to introduce extra 
terms into the system density matrix. As the system gets 
near an ionization threshold, the system becomes partic- 
ularly delicate, since the particle with the highest energy 
can tunnel further away from the center of the trap and 
thus external disturbance has bigger exciting effect on 
the system. Moreover, immediately after the ionization 
threshold is passed, the system density matrix is sub- 
ject to various excitations due to wavefunction collapses. 
These excitations are crucial to the fidelity of Fock state 
production, since they cause significant reversion in the 
number uncertainty of the final state. A characteristic 
measurement of tendency of excitation is the energy gap, 
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FIG. 5: (color online) Excitation energy gaps between ground 
and first-excited states as function of trap depth for 2 (circle) , 
3 (square), 4 (diamond), 5 (upright triangle), and 6 (invert 
triangle) sodium atoms. Trap size is 5 rim. 
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A, which is define as the difference between total energies 
of ground and first-excited Bethe ansatz states (if both 
exist). According to our calculation, they are on the or- 
der of a few ks X 10nK (see Fig. [5]). A puts restrictions 
in two-fold. Firstly, the temperature must be maintained 
lower than a few lOnK, otherwise, fidelity could be en- 
dangered due to thermal excitation. Furthermore, the 
energy gap puts a requirement on the adiabaticity con- 
dition [24j : the culling speed must be much smaller than 

A 2 

=jr- to maintain a relatively high fidelity. To give a sense 
of number, we consider the culling of trapping potential 
from the ionization threshold of 3 particles down to that 
of 2 particles at trap size of 5/xm and transverse trap- 
ping frequency lj±_ = 2n x 150kHz. According to our 
calculation, the minimum time required to complete this 
portion of the culling should be no less than 0.3ms to be 
considered as adiabatic. 



V. CONCLUSION 

In conclusion, we calculated the conditions for num- 
ber states of ultracold atoms in ID optical trap with the 
Bethe ansatz approach. We charted ionization thresholds 
in the parameter space. We also discussed the quantum 
mechanical processes in producing number states in the 
ideal case and the effect of excitations. 
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APPENDIX A: VALID BETHE ANSATZ 
SOLUTIONS 

First of all, we need a change of unit to make the de- 
pendence of secular equation on interaction strength ex- 
plicit. (See Section U for the previous choice of unit.) To 
that end, we choose the trap length L as the length unit, 
H 2 /mL 2 as the energy unit. Then the secular equation 
Eq. (fTT)) is transformed to 

-^1 = ^+2 sin" 1 + 

g [ tan -i (k±h^j + tan -i , ( ai) 

Now Eq. (|A1|) depends both on the set of integers I and 
the interaction strength c. For the mere purpose of solv- 
ing Eq[17l any set of integers / may be used. Note that 



the choice of integers is discrete and is usually enumer- 
able while that of the interaction strength c is continu- 
ous and non-enumerable. Based on this fact, we claim 
that if a given set of integers lead to valid solution at 
some interaction strength, it does so at any interaction 
strength. Particularly, the solution in the weak interac- 
tion region should approach that of the non-interacting 
case as c — > 0. Solution of the latter is a well-taught 
exercise in many quantum mechanics textbooks. 

Thus, a 'promising solution' to Eq. ljAlj) should con- 
verge onto that of the c — case and we can use the 
known solutions to reject spurious 'solutions' for the in- 
teracting cases. With a few examples of N and upper 
limits M , we exhaust all the combinadics of N numbers 
from the range [0, M] and solve Eq. (|Al[) with each com- 
bination. Our experiments show that in the limit c — > 0, 
a solution approaches that of c = if and only if the 
set / consists of positive and mutually distinct integers. 
Furthermore, we also find that the wave numbers in the 
solution are mutually distinct if and only if the integers 
in the set / are mutually distinct. 



APPENDIX B: ORDER OF BETHE ANSATZ 
STATES 

With the valid solutions found in Appendix [A] it still 
left to determine which one is the ground state and which 
is the first excited state, and so on. In a similar manner 
as in Appendix [AJ we will argue based on intuition that 
for any given N, the set I = {1, 2, • ■ • , N} leads to the 
ground state and the set / = {1, 2, • ■ • , N — 1, N + 1} to 
the first excited state. 

Our clue comes from the strong interaction limit. As 
is well known, in the strongly interacting limit the par- 
ticles behave as fermions [9(. Thus the ground state of 
our system must be like that of a degenerate fermion 
system. Our numerical calculations show that solving 
Eq.flATJ with the set / = {1,2,- ■■ , N} leads to a so- 
lution with energy that approaches that of the ground 
state of the degenerate fermion system in the limit c — > 

00. Furthermore, the solution obtained with the set 
I = {1, 2, • • ■ , N — 1, N + 1} approaches that of the first- 
excited state of the same system in the same limit. We 
thus established what are ground state and first excited 
state for our system in general. However, there is little 
we could say beyond that. Within limit of calculation er- 
ror, our experiment is not conclusive about which is the 
second excited states. It might be 1, 2, • • • , N — 1, N + 2, 

1, 2, • • • , N - 2, TV, N + 1, or still other, depending on N 
and the trap parameters. In general the order in the en- 
ergy level of a Bethe ansatz state depends on both the 
highest quantum number and the total of these quantum 
numbers. For complete ordering, one need something as 
what Hund's Rule is in atomic physics. For our paper, 
it suffices to know just the ground state and first excited 
state. 
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